Two dimensional Symmetry Protected Topological Phases with PSU(iV) 

and time reversal symmetry 
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Symmetry protected topological phase is one type of nontrivial quantum disordered many-body 
state of matter. In this work we study one class of symmetry protected topological phases in two 
dimensional space, with both PSU(iV) and time reversal symmetry. These states can be described 
by a principal chiral model with a topological 0— term. As long as the time-reversal symmetry and 
PSU(iV) symmetry are both preserved, the 1+1 dimensional boundary of this system must be either 
gap less or degenerate. We will also construct a wave function of a spin-1 system on the honeycomb 
lattice, which is a candidate for the symmetry protected topological phase with both SO(3) and 
time-reversal symmetry. 

PACS numbers: 



1. INTRODUCTION 

The interplay between strong interaction and strong 
quantum fluctuation can lead to many remarkable prop- 
erties in quantum disordered phases. In the most trivial 
case, a quantum disordered phase is fully gapped, and 
nondegenerate, and its ground state wave function can 
be adiabatically connected to a direct product wave func- 
tion without any phase transition. These quantum dis- 
ordered phases are "trivial" , in the sense that they are 
quantum analogues of classical disordered phases, namely 
they are completely featureless. The best example of triv- 
ial quantum disordered phase is the Mott insulator phase 
of spinless bosonic atoms trapped in an optical lattice. 
Every state of this Mott insulator can be adiabatically 
connected to direct product state Yii \ fk — tyu where 
hi is the boson number operator on site i, and k is an 
integer. 

The description of trivial quantum disordered phases 
is semiclassical, i.e. we can describe the trivial quantum 
disordered phase using a field theory defined with a Lan- 
dau order parameter. For example, a trivial quantum 
disordered phase of a system with SO(3) spin rotation 
symmetry can be described by a semiclassical nonlinear 
sigma model (NLSM) defined with the Neel order vector 
n with unit length (n) 2 = 1: 



S = J d d xdr -(d^n) 



(1) 



In spatial dimensions higher than 1, by tuning the pa- 
rameter g, there is an order-disordered phase transition: 
when g < g c , n is ordered, and the SO(3) symmetry 
is spontaneously broken; when g > g c , n is disordered, 
and the ground state wave function of this quantum dis- 
ordered phase is approximately a direct product state 
0) ~ Yii 1^ = where I is the angular momentum 
quantum number on every site. 

Now there is a consensus that, quantum disordered 



phases of strongly correlated systems can have much 
richer structures compared with classical disordered 
phases. Roughly speaking, three types of "nontrivial" 
quantum disordered phases have been studied: (i.) topo- 
logical phases, whose ground state is fully gapped but 
topologically degenerate; (ii.) Algebraic spin (Bose) liq- 
uid phase, which is still a quantum disordered phase of a 
bosonic spin system, but the spectrum remains gapless, 
and the physical quantities have a power-law instead of 
short range correlation; (in.) Symmetry protected topo- 
logical (SPT) phases, whose bulk spectrum is identical 
to a trivial quantum disordered phase, but its boundary 
must be either gapless or degenerate when and only when 
the system has certain symmetry G. The ground state 
wave function of a SPT is completely different and can- 
not be continuously connected to a trivial wave function 
without a phase transition, when the Hamiltonian is in- 
variant under symmetry G. The 2d quantum spin Hall 
insulator and the 3d topological insulator are both SPT 
phases with time-reversal symmetry. 

In this work we will focus on SPT phases of bosonic 
spin systems. SPT is a pure quantum phenomenon, and 
there is no analogue in classical world. Thus one would 
expect that the description for these states should be 
purely quantum, and a semiclassical formalism should 
completely fail to describe them. However, we will try 
to demonstrate that, the SPT phases can still be de- 
scribed semiclassically using NLSMs like Eq. [1] as long 
as we include an appropriate topological term. We want 
to stress that, in our approach, the target space of the 
NLSM is the manifold of a semiclassical orderparameter. 
This is different from the approach in Ref. 0, Q , where 
a NLSM was also introduced to describe SPT phases, 
but the target space of this NLSM is the group manifold 
of the symmetry G, instead of the manifold of an order 
parameter. 

At least in one dimensional systems, semiclassical 
NLSMs have been proved successful in describing SPT 
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phases. For example, the 0(3) NLSM in Eq. [TJ plus a 
topological O— term describes a spin-1 Heisenberg chain 
when = 27r, and it is well-known that the spin-1 an- 
tiferromagnetic Heisenberg model is a SPT phase with 
2-fold degeneracy at each boundary 0, Q ■ This two fold 
degeneracy at the boundary can be read off from this Id 
0(3) NLSM, since its boundary is a O+ld 0(3) NLSM 
with a Wess-Zumino-Witten term at level k = 1, whose 
ground state is two fold degenerate 

In principle, a NLSM describes a system with long 
correlation length. Thus a NLSM plus a 0— term most 
precisely describes a SPT phase tuned close to a crit- 
ical point (but still in the SPT phase). When a SPT 
phase is tuned close to a critical point, the NLSM not 
only describes its topological properties (e.g. edge states 
etc.), but also describes its dynamics, for example excita- 
tion spectrum above the energy gap (much smaller than 
the ultraviolet cut-off). When the SPT phase is tuned 
deep inside the SPT phase, namely the correlation length 
is comparable with the lattice constant, this NLSM can 
no longer describe its dynamics accurately, but since the 
topological properties of this SPT phase is unchanged 
while tuning, these topological properties can still be de- 
scribed by the NLSM. 

In Ref. [6], the author discussed a class of three di- 
mensional SPT phases with SU(iV) symmetry (Rigor- 
ously speaking, the symmetry of the SPT constructed in 
Ref. @ is PSU(A) = SV(N)/Z N , where Z N is the center 
of SU(A)), and just like the Haldane phase in Id, these 
3d SPT phases are described by a NLSM defined with the 
Neel order parameter only. For SU(7V) magnet, the man- 
ifold M of the Neel order is M = u(m) ^_ m) @> H- 
When N > 2 and < m < N, tt 4 [M] = Z. Thus a non- 
trivial topological 0— term can be defined for the SU(7V) 
Neel order parameter, and when = 2tt, it was argued 
that the system is a three dimensional SPT, whose 2+ld 
boundary must be either gapless or degenerate jg]. 



In 2 dimensional space, a straightforward general- 
ization of Ref. [f| is difficult, since for N > 2, 

7r 3[ U ( m )xU('iv- m ) ] = Zl ' tnus a e ^ term is not wel1 de- 
fined for manifold M. in two dimensions with general N. 
For the simplest case with N = 2 and m = 1, i.e. the or- 
dinary SU(2) Neel order whose manifold is S 2 , although 
the homotopy group Tts[S ] = Z, a nontrivial mapping 
from the space-time manifold to the target space S 2 can- 
not be written as an integral of local terms of the Neel 
order parameter, thus it is much more complicated. Al- 
ternatively, in this paper we will study 2d systems with 
symmetry PSU(iV) x Zj, where is the time-reversal 
symmetry. We will demonstrate that for a system with 
PSU(iV) x Z2 symmetry, a topological 0— term can be 
defined in 2+1 dimensional space-time with semiclassical 
order parameters, and when = 2tt the topological term 
will drive the system into a SPT phase. 



2. FIELD THEORY DESCRIPTION 

2.1 SPT phase at 9 = 2tt 

In this section we will discuss the field theory descrip- 
tion of the 2+ld SPT phase with PSU(iV) x Z^ sym- 
metry. As we discussed in the introduction, we will take 
the semiclassical formalism, and define the field theory 
in terms of the order parameter of PSU(N), whose con- 
figurations form manifold M = U(m ) x |^_ m ) ■ Every 
element V of the manifold M. can be represented as 

v = v^nv, n = I 1 , (2) 

OjV— mxro> liV— mxN— m 

where V is a SU(7V) transformation matrix. V is a Her- 
mitian matrix that satisfies V 2 = Inxn- In fact, when 
N = 2, M. is precisely S 2 . The matrix order parameter 

V is always invariant under the center Zm of SU(7V): 

V = e i2 * k / N l NxN (k = 1,---N - 1), thus a NLSM 
defined with V has symmetry PSU(A^) = S\J(N)/Z N . 
When N = 2, PSU(2)=SO(3) is the ordinary spin rota- 
tion group of model Eq. [TJ If the symmetry of a quantum 
spin system is SO (3) instead of SU(2), then the Hilbert 
space on every site of the system must be a representation 
of SO (3), thus we are restricted to integer spin systems 
only. With integer spins, (Zj) 2 = +1. 

For a general N, The homotopy group ^[A^] = Zi, 
thus a NLSM defined with V does not have a nontriv- 
ial topological 0— term in 2+ld space-time. Thus for a 
2+ld system with PSU(iV) symmetry only, there is no 
nontrivial SPT that can be described using semiclassical 
order parameter V. Now let us combine PSU(TV) and 
time-reversal symmetry together, and define the follow- 
ing order parameter U : 



U = cos(6»)1jvxw + ism{9)V. 



(3) 



Now U is a unitary matrix, and U G U(iV). Since 
7T3[U(7V)] = Z, a principal chiral nonlinear sigma model 
can be defined with order parameter U, plus a nontrivial 
topological 0— term: 



S = J d 2 xdr ^ixldnrfdpU] 
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tr [C/ f <9 M UU^ d v UU^ d p U] e^p 



(4) 



In the simplest case with N = 2 and m = 1, the order 
parameter V is equivalent to an 0(3) vector ft: V = n-B. 
Then Eq.His equivalent to an 0(4) NLSM: 



(5) 



S = I d 2 xdr -(dp 



ie 
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where = (cos(#), s'm(6)n). 

The matrix U has a SU(A)i e f t transformation and a 
SU(A) right transformation, but there are higher order 
terms in the Eq.[l]that break the two SU(iV) symmetries 
down to its diagonal subgroup PSU(iV). For general N, 
we can define the following transformations: 



SU(A) : P^V^PV, 



V -> V*, 



->• 7T — 



V^V*. (6) 



Under this definition, the SU(iV) and transformations 
commute with each other. 

In the follows, we will argue that, when 9 = 2tt, Eq. [4] 
with symmetry SU(iV) x describes a SPT, whose 1+ld 
boundary must be either gapless or degenerate. 

In Eq. |4j by tuning g, there is an order-disorder phase 
transition. We will always focus on the disordered phase 
of Eq. 21 thus we will focus on the phase with a large 
coupling constant g. When 9 = 2tt, the bulk spectrum 
of the quantum disordered phase is identical to the case 
with 9 = 0, thus the bulk is fully gapped and nondegen- 
erate. Then one can safely integrate out the bulk and 
focus on the boundary. At 9 = 2tt, the boundary of 
the system is described by the following principal chiral 
model with a Wess-Zumino-Witten (WZW) term: 



S b = J dxdr \x[d^U ] d,JJ] 



dxdrdu 



i2n 
24tt 2 



Here U(x, r, u) is an extension of physical order parame- 
ter U (x, r) that satisfies 



U(x, r, u = 0) = 1 



NxN, 



U{x,r,u = 1) = U{x,r). 



(8) 



For the simple case with N = 2m = 2, Eq. [7] reduces to 
a 1+ld 0(4) NLSM with a WZW term at level- 1: 



S 



dxdr -(d^Y 



dxdrdu 



2-ki 
12^2 



iabcd^ p rd^ b d^ c d p 4 d , (9) 



This principal chiral model Eq. with a full 
SU(iV)i e ftxSU(-/V) r ight symmetry, was proved to be a 
gapless conformal field theory [{J, However, in our 
system the SU(A)i G f t xSU(A) r ight symmetry is broken 
down to the diagonal SU(iV), thus the conformal field 
theory might be gapped out due to relevant perturba- 
tions introduced by this symmetry reduction. With this 
symmetry reduction, the boundary theory Eq. [7J is re- 
duced to the following NLSM with a 9'— term: 

S b = [ dxdr ltr[(d^T) 2 } + ^tr[Vd^Vd v V]e^, (10) 

J g io7T 



as long as the Z^ symmetry is preserved, namely the 
expectation value of cos(0) is zero, the boundary 9' 
is fixed at 9' = tt. Under the Z|" transformation, 
9' — > 2tt — 9'. If the time- reversal symmetry is explicitly 
broken, namely a background field that couples linearly 
to cos(#) is turned on, then at the boundary 9' is also 
tuned away from tt: 9' = 29 — 2cos(6>) sin(6>). 

If we ignore the physical interpretation of the field P, 
this 1+ld NLSM at 9' = tt (Eq. HHJ) can be used to 
describe the SU(iV) antiferromagnet with conjugate rep- 
resentations on A and B sublattices 0, @, Ef , and V is 
precisely the SU(7V) Neel order parameter. In fact, for 
the simplest case with N = 2, m = 1, Eq. [TOl precisely 
reduces to an 0(3) NLSM with a 9' term with 9' = n: 

S = ( dxdr -(d^ri) 2 + l -^-e a bc£^n a d f _ l n b d„n c . (11) 

J 9 87T 

With 9' = tt, this model describes an antiferromagnetic 
spin-1/2 chain, and based on the Lieb-Schultz-Mattis 
(LSM) theorem this theory must be either gapless or de- 
generate [l2| . 

In Eq. 1101 9' = n is the transition point between two 
stable fixed points at 9' = and 2tt, and this transi- 
tion is driven by tuning 9'. For example, when m = 1, 
Eq. [TU1 becomes the CP^ -1 model with 9' = tt, and in 
the large— A'" limit it was demonstrated that the CP^ -1 
model has a first order transition at 9' = tt ll|, Il3 |. 



For general m and N, this transition at 9' = tt was 
analyzed through a renormalization group calculation of 
both g and 9' as in Ref. [IMS- If this transition is 
continuous, then this 1+ld boundary system must be a 
gapless CFT at 9' = tt; if this transition is first order, 
then this boundary system must be two fold degenerate 
at 9' = tt 19]. Thus we conclude when the bulk theory 
Eq. @]has 9 = 2tt, its boundary must be nontrivial, i.e. 
it must be either gapless or degenerate. 

To demonstrate that this phase is a SPT, we need 
to show that its boundary can only be realized as the 
boundary of a 2+ld system, i.e. it cannot be realized as 
a real one dimensional lattice quantum spin system with 
the same symmetry. For example, the boundary of a 2d 
quantum spin Hall insulator is a Id helical Luttinger liq- 
uid with central charge 1, and it can be argued that this 
Id helical Luttinger liquid cannot be realized as a Id elec- 
tron system with time-reversal symmetry (20j . Also, the 
boundary of 3d topological insulator is a single 2d Dirac 
cone, which cannot be realized in a pure 2d system with 
time-reversal symmetry. Let us take the simplest case 
with N = 2m = 2 as an example. In order to argue that 
Eq.rTTl cannot be realized in a Id system with SO(3)xZ^ 
symmetry, let us break the time-reversal symmetry at the 
boundary, but make a domain wall of the time-reversal 
symmetry breaking pattern: 



0' > tt, for x > 0; 9' < tt, for x < 0. 



(12) 
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The two sides of the domain wall are conjugate to each 
other under time-reversal transformation. Based on the 
renormalization group calculation of Ref. 14-lj|, and 
the argument in Ref. when 6' > 7r the system is 
in the same phase as G' = 2tt, while when ©' < tt the 
system is in the same phase as 6' = 0, both phases are 
fully gapped and nondegenerate. At the domain wall this 
system is described by a 0+ld 0(3) NLSM with a WZW 
term at level- 1: 

/l f 2ni 

dr -{d^nf + J drdu —t tiv t abc n a d fl n h d v n\l'i) 

which is precisely the model describing a single spin- 1/2. 
A spin-1/2 excitation is not a representation of SO(3) 
(it is a representation of SU(2)), since it is not invari- 
ant under the center of SU(2). Also, under the square 
of time-reversal transformation, the wave function of a 
spin-1/2 excitation would change sign: (Z^) 2 = —1. This 
implies that physical symmetries fractionalize at the do- 
main wall. 

Although one dimensional spin chains can have frac- 
tionalized excitations, this phenomenon of deconfincd 
domain- wall fractionalization cannot happen in a one di- 
mensional integer spin chain. Consider for example two 
Id systems with integer spins and SO(3) symmetry, then 
if these two systems are conjugate to each other under 
time-reversal symmetry, then they must be either both 
Id SPT, or both trivial states. Then at their domain wall 
there should be either an integer localized domain wall 
spin, or no domain wall spin at all. 

The analysis of domain wall state at the boundary can 
be generalized to arbitrary N. For arbitrary N , the the- 
ory Eq. [4] is parametrized by integer m. For general N 
and m, a single spin with representation Fig. [IJ> is lo- 
calized at the domain wall Eq. [121 an d the theories with 
different m (mod N) have inequivalent domain wall spins. 
This domain wall state is always a fractionalized excita- 
tion of PSU(TV), since it is not invariant under the center 
of SU(iV). Based on all the analysis above, we can con- 
clude that Eq. [4] with = 2tt is a SPT with symmetry 
PSU(7V)xZ!f. Our formalism suggests that for general 
N, there are at least N inequivalent phases: there is one 
trivial phase, and N — 1 SPT phases described by Eq. [J 
with m = 1, • • - N— 1, which have different localized spins 
at the domain wall Ref. [12] For the simplest case with 
N = 2m = 2, i.e. the case with SO(3)xZ^ symmetry, 
there are only two inequivalent phases, which is consis- 
tent with the classification in Ref. [l|. For general N, our 
prediction can be compared with future group cohomol- 
ogy computation. 



2.2 Physics with 6 = tt 



exp(i"0/2) 




i] 



B 



FIG. 1: (a) We map Eq. \T7\ to a one dimensional tight- 
binding model. Hopping between nearest neighbor sites cor- 
responds to changing $ by 1/2. The 0-term grants two types 
of monopoles a factor exp(iO/2) and exp(— iQ/2) respec- 
tively, which forbids nearest neighbor hopping when O = n 
due to destructive interference between these two types of 
monopoles. (6) The SU(7V) spin localized at the domain wall 
Eq.ETJ] 

generacy. In Ref. 0, using a different argument, it was 
concluded for general principal chiral models Eq. [4] that 
their quantum disordered phases must be either gapless 
or two fold degenerate when O = tt. In Ref. [19] it was 
assumed that the system has a full SU(iV)i e ft x SU(7V) r i g ht 
symmetry. In our current case, the actual symmetry is 
PSU(iV)xZj. In this section, using a different argument 
from Ref. [l9j, we will make the same conclusion for Eq.[H 
i.e. its quantum disordered phase cannot be gapped with- 
out degeneracy when O = tt. 

In order to argue that a system is either gapless or de- 
generate when O = 7r, we only need to argue that if the 
system is gapped, it must be degenerate. Thus let us as- 
sume it is gapped in the disordered phase of O = tt. Un- 
der this assumption, the coupling constant g must flow to 
infinity in the infrared limit under renormalization group, 
this is because if g flows to any fixed point with finite con- 
stant g* , then the system must be scaling invariant and 
gapless. Thus g must flow to infinity once we assume the 
system is gapped. Now let us take g to infinity, then the 
first term of Eq. [4] vanishes, and Eq. [4] reduces to a pure 
topological 6— term: 



S = I d xdr 



ie 
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txltfd^UtfduUtfdpU^up. (14) 



It was first discussed in Ref. [2l| that, the 0(4) NLSM 
Eq. [5] at O = 7r cannot be trivially gapped without de- 



This O— term contributes phase factor exp(iO) to every 
SU(7V) instanton in the space-time. However, since our 
system only has PSU(A) x Z^ symmetry instead of a full 
SU(7V) xSU(^V) symmetry, an instanton will fractionalize 
into two monopoles. A monopole centered at the origin 
has the following configuration: 

U = cos(9)l NxN + isin(9)P, 

9(\R\ = 0) = 0, or tt, 9(\R\ = oo) = tt/2; 
f <? R J— tT [ Vd VduV ] e = i. (15) 

J\R\=R 167r 

Here R is the coordinate in the 2+ld Euclidean space- 
time. In the simplest case with N — 2m — 2, this 
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monopole is the ordinary "hedgehog" monopole of vec- 
tor ri in the space-time. When 9 = and tt at the origin 
\R\ = 0, this monopole carries instanton number 1/2 and 
— 1/2 respectively, thus this monopole contributes phase 
factor exp(±i0/2) to the partition function. 

Since we are interested in the bulk physics, we can 
compactify the two dimensional space into a sphere S 2 . 
Now let us define the following quantity $(r) for every 
time slice r: 

$(r) = f d 2 x -^trlVdiVdjVfoj, (16) 
J oZir 

and since ^[.A/l] = Z, <E>(t) is quantized as integer or 
half-integer, as long as the configuration of V has no 
singularity at time r. $(t) is increased and decreased 
through the monopoles described in the previous para- 
graph, and one monopole in the space-time changes $ by 
1/2: $(r = +00) — $(— 00) = n m /2, where n m is the 
monopole number in the 2+ld space-time. 

Now under the assumption that the system is gapped 
(hence g flows to infinity) , Eq. H] and Eq. Q3] reduce to 
the following single particle quantum mechanics problem 
defined on a periodic lattice with lattice constant 1/2 
(Fig. Eh): 

S= J dr im(<9 T $) 2 ±iOd T $ + V(<f>), m ->• 0, (17) 

where V(&) is a deep periodic potential that makes <£> 
takes only integer and half-integer values. Thus the orig- 
inal principal chiral model has now reduced to a single 
particle tight binding model, where each lattice site corre- 
sponds to a quantized value of Hopping from site q to 
site q+ 1/2 corresponds to a monopole in the space-time, 
and there are two different types of monopole, depending 
on the sign of cos(9) at the monopole core. The 9-term 
will contribute a factor exp(iO/2) and exp(— iO/2) to the 
monopole with cos(6>) = +1 and cos(#) = — 1 at the core 
respectively. With time-reversal symmetry, the two types 
of monopoles are degenerate, thus when O = tt these two 
types of monopoles have complete destructive interfer- 
ence with each other, i.e. hopping by one lattice constant 
in Fig. QJi is forbidden. However, hopping by two lattice 
constants is still allowed, but the band structure will be 
doubly degenerate, namely on this one dimensional lat- 
tice (Fig. [T}i) the states with lattice momentum p = 
and p = 2ir are degenerate. 

If is tuned away from tt, then the nearest neighbor 
hopping in the tight-binding model Eq. [17] is allowed, 
and the ground state of Eq. [17] is nondegenerate. Now 
we have argued that once the system is gapped at = tt, 
it must be two fold degenerate, namely the system must 
be either gapless or degenerate at = tt. The analysis in 
this section implies that when we tune = 2tt to 0, the 
bulk spectrum must change at = tt, i.e. the SPT phase 
and trivial phase must be separated by a bulk transition 
at = tt. 




spin-l/2 



FIG. 2: Lattice construction for N = 2, m = 1 and its edge 
state, (a) the mean field Hamiltonian Eq. consists of the 
nearest-neighbor hopping t and the "color" -orbit interaction 
A and t'. The dashed lines represent the usual "color" -orbit 
coupling ~ A in the Kane-Mele model. It is equivalent to 
the Kane-Mele model except for an additional x-directional 
"color" -orbit couplings ~ t' (represented by the dotted lines) 
on top of the usual "color"-orbit term ~ A. This anisotropic 
coupling breaks the gauge symmetry down to SU(2) and the 
lattice rotational symmetry which is irrelevant for the sym- 
metry protected phase, (b) the edge theory consists of right- 
moving (t z = 1) and left-moving (r z = —1) spin-l/2 dou- 
blets. The time-reversal symmetry maps a right-mover into 
a left-mover with the opposite spin state, i.e. the color and 
the spin are both flipped by the time-reversal symmetry, (c) 
the domain wall in <j>o(x) Eq.(24) traps a localized spin-l/2. 
This can be derived from the WZW term written in terms of 
(4>o, 0/j), /i = 1, 2, 3 as discussed below Eq. 1111 



3. LATTICE CONSTRUCTION FOR SPT WITH 
SO(3)xZ^ SYMMETRY 

In this section we will try to construct a lattice spin 
state for the 2d SPT with SO (3) and time-reversal sym- 
metry. Since the Hilbert space on every site must be a 
representation of SO (3) group, we must start with an in- 
teger spin system. Let us consider a spin-1 system on a 
honeycomb lattice, and we will construct a spin many- 
body wave function using the following two-color slave 
fcrmion formalism, which was introduced to understand 
the spin liquid phenomena observed in Ba3NiSb20g [22| : 

= 2 S ^2 f»-a.,i a apfl3,a,i- ( 18 ) 
a,/3=t,4- o=X,2 

Here cr^ are three spin-l/2 Pauli matrices. Each spinon 
/ Q , Q has two indices: a =T, 4- denotes spin and a = 1, 2 
is a "color" quantum number. Thus we can consider not 
only the usual spin SU(2) rotations in the a — /3 space, 
but also color SU(2) transformations in the a — b space. 
Matching with the spin Hilbert space requires not only 
constraining the total fermion number to half-filling (two 
fermions per site), but also requiring each site to be an 
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color SU(2) singlet, which guarantees that on each site 
the spin space is a symmetric spin-1 representation: 

hi = E E fa,a,if a > a >i = 2, 
o=l,2a=ti 

f M = E fla,i< b fa, b ,i = 0. (19) 



Written in terms of the Majorana fermion ry, the mean 
field Hamiltonian reads: 

H M f = ^2 -tr)lTx 2 rij + ^ iXi'. j/'l 

<i,j> «,j> 

+ ^ v:), • // ''-- ( 25 ) 



Here r£ b are three Pauli matrices that operate on the 
color indices. Under time-reversal transformation, in or- 
der to satisfy the commutation relation between Pauli 
matrices, r M should transform just like spin operators: 

Due to these two independent constraints in Eq. 1191 
the spinon f a a will have a gauge symmetry. In order 
to identify the full gauge symmetry, we need to rewrite 
fa,a,i in terms of Majorana fermions rj as follows: 



(20) 



On every site, r\i has in total three two-component spaces, 
making the maximal possible transformation on r\i SO(8). 
Within this SO(8), the spin SU(2) transformations are 
generated by the three operators 



(a x X y , a y , a z X y ), 



(21) 



where the Pauli matrices A a operate on the two- 
component space (Re[/], Im[/]). Under time-reversal 
transformation, r\ — > a y T y X z rj 1 i.e. both spin and colors 
are flipped under Zj. 

The total gauge symmetry on rj is the maximal sub- 
group of SO(8) that commutes with the spin-SU(2) op- 
erators. This is Sp(4) ~ SO(5) generated by the ten ma- 
trices r a6 = jj[r a ,r&], where 



- a y T y x z , r 2 = a y r y x x , r 3 = t v x v , 



r. 



(22) 



These T a with a = 1 • • • 5 define five gamma matrices that 
satisfy the Clifford algebra {r a ,Ph} = 28 a b- Tah and T a 
are all 8 x 8 hermitian matrices. r a & are all antisymmetric 
and imaginary, while T a are symmetric. 

A spin wave function can be constructed through a 
slave fermion wave function, by projecting the mean field 
ground state to satisfy the gauge constraints: 



|G spi n) =n p ("* = 2 )®P(ff =0)|/ Q , a > 



(23) 



Now let us consider the following mean field Hamilto- 
nian of slave fermion on a honeycomb lattice: 



H MF = Y. ~*/</i+ X iXu vfla T abfj,b 

<i,j> <i,J> 



When t' = 0, this mean field Hamiltonian is a quan- 
tum "color" Hall Hamiltonian, i.e. it is equivalent to the 
Kane-Mele quantum spin Hall Hamiltonian |23l . 24 1, al- 
though instead of a spin-orbit coupling, in Eq. [24] the A 
term is a color-orbit coupling. At the mean field level, the 
quantum color Hall mean field Hamiltonian Eq. [24] has 
edge states: there is a left-moving spin- 1/2 doublet slave 
fermion with t z = 1, and another right-moving spin-1/2 
doublet slave fermion with r z = — 1. Backscattering be- 
tween left and right moving slave fermions is forbidden, 
as long as the time-reversal symmetry and spin rotation 
symmetry are preserved. 

The t' term is a color-orbit coupling between 2nd 
neighbor sites along the x directions only. If t' = 0, 
the color-orbit coupling term A breaks the Sp(4) gauge 
symmetry down to its subgroup SO (4), which is gener- 
ated by T a b, with a, b = 1, • • • 4; when t' is nonzero, the 
gauge symmetry is broken down to SU(2) generated by 
T a b, with a, b = 1, 2, 3. Since the t' term is time-reversal 
even, when t' is small compared with t and A, the edge 
state cannot be gapped out without degeneracy. Now the 
edge states can be describe by the following 1+ld field 
theory: 



3 



(26) 



where A 1 ^ is the residual SU(2) gauge field, and G a — 
e a b c Tbc: a,b,c= 1, 2, 3. r/ is the Majorana fermion intro- 
duced in Eq. 7 = r y , 7 1 = r x , 7 s = t 2 , fj = 77*7°. 

Eq. [26] is precisely the same field theory that describes 
the spin-1/2 chain, if we take the standard SU(2) gauge 
field formalism for spin-1/2 systems 25[. There is no 
symmetry allowed fermion bilinear terms that can be 
added to Eq. [26] It is well-known that a spin-1/2 chain 
must be either gapless or degenerate, and when it is gap- 
less it can be described by the 1+ld SU(2)i conformal 
field theory (plus marginal perturbations). Thus we con- 
clude that the boundary of the bulk state Eq. [M] is ei- 
ther a gapless SU(2)i CFT, or degenerate due to spon- 
taneous time-reversal symmetry breaking. The sponta- 
neous time-reversal symmetry breaking can be induced 
by a relevant four fermion interaction term in Eq. 1261 

It is well-known that the SU(2)i CFT (and spin-1/2 
chain) is equivalent to an 0(4) NLSM with a WZW term 
(Eq. gjj. The WZW term of Eq. can be derived from 
Eq. [55] by coupling the fermions 77 in Eq. [5S] to a four 
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component order parameter <j> with unit-length: 



Doll 



(27) 



k=l 



where S k are the three matrices defined in Eq. [21] which 
generate the spin rotations. Careful analysis shows that 
the order parameter <fr defined here has the same transfor- 
mation as the 0(4) vector in Eq.[SJ With the coupling in 
Eq. [27J after integrating out the slave fermions, a WZW 
term at precisely level- 1 (Eq. [9]) will be generated [2r| . 
Thus the field theory discussed in the previous section 
can be precisely derived from this lattice construction. 
<po changes sign under time-reversal transformation. If 
a domain wall of 4>o is created, a spin- 1/2 excitation is 
localized at the domain wall, consistent with the field 
theory discussed below Eq. [TT] (see also Fig. [2]). 

In the bulk the slave fermion is gapped. Since the time- 
reversal symmetry in the bulk excludes the existence of 
a Chern-Simons term for the residual SU(2) gauge field, 
this nonabelian gauge field will lead to confinement, and 
this confined state has no topological degeneracy. At the 
boundary, the effect of this confinement is not totally un- 
derstood. For instance, this confinement might gap out 
the boundary state through a spontaneous time-reversal 
symmetry breaking, namely the order (7777) ^ is sponta- 
neously generated. But nevertheless, the boundary will 
not be gapped out without degeneracy. 

This lattice construction of SPT can be checked nu- 
merically in the future. Let us define the system on a 
torus without boundary, and a spin wave function can 
be constructed by gauge projecting the mean field state 
Eq. 1241 Given this projected spin wave function, one can 
numerically compute various quantities such as spin-spin 
correlation, topological entanglement entropy, and entan- 
glement spectrum. Since the bulk is completely gapped, 
the spin correlation should be short ranged. And since 
the bulk of the system has no topological degeneracy, the 
topological entanglement entropy, which is defined as a 
universal constant term in addition to the standard area 
law, should be zero. However, since the system has non- 
trivial edge states, the edge states should also exist in 
the entanglement spectrum, which can be computed us- 
ing the projected wave function. We will leave these to 
future studies. 



The principal chiral model Eq. [4] can describe some 
other symmetry protected topological phases as well. For 
example, the spin-2 AKLT phase on the square lattice is 
known to have a nontrivial Id edge states. Unlike the 
Id Haldane phase, AKLT states at higher dimensions 
require translation symmetry to protect its edge states. 
For example, on the square lattice, the edge state of the 
spin-2 AKLT phase is a spin-1/2 chain, and if the trans- 
lation symmetry is broken, this edge spin-1/2 chain will 
be dimerized and gapped. The spin-2 AKLT phase on 
the square lattice can be viewed as a SPT phase with 
SO(3)xZ2 symmetry, while here the Z2 is translation by 
one lattice constant instead of time-reversal transforma- 
tion. Then Eq.[5]and Eq.[4]with 6 = 2ir can also describe 
the two dimensional AKLT phase, and its generalizations 
to arbitrary N. 

Our result and analysis apply for all even spatial di- 
mensions. For large enough N and m, 7r2d+i[SU(iV)] = 

Z, and 7r2d+i[ u(m) xu(w-m) ] = Zl " Thus at least based 
on the field theory, a SPT with PSU(JV)xZf symmetry 
exists in arbitrary even spatial dimension, and it is al- 
ways described by a principal chiral model defined with 
order parameter U introduced in Eq. [3J 

Besides a general classification given in Ref. I], SPT 
phases have attracted a lot of attentions and efforts re- 
cently 27H31j . Most of these studies are focused on SPT 
phases with a U(l) symmetry that is associated with bo- 
son number conservation, and with a U(l) symmetry it 
is convenient to describe the system using a dual U(l) 
gauge field in 2+ld, and a topological response function 
of electromagnetic field in 3+ld. The SPT phases dis- 
cussed in the current paper and in our previous study 
(Ref. Q) may also have a description in terms of their 
responses to external background gauge fields. We will 
leave this to future study. 
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4. DISCUSSION AND SUMMARY 



In this work we studied a class of two dimensional sym- 
metry protected topological phases with PSU(7V)xZj 
symmetry. These SPT phases are described by a 2+1 
dimensional principal chiral model with O = 2ir (Eq. [4]) , 
and their boundary states are described by a 1+ld NLSM 
with ©' = tt, which must be either gapless or degenerate 
when the symmetry PSU(-/V)xZ^ is preserved. 
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